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Abstract 



A QCD analysis is reported of ATLAS data on inclusive W^ and Z boson production in pp collisions 
at the LHC, jointly with ep deep inelastic scattering data from HERA. The ATLAS data exhibit sensitivity 
to the light quark sea composition and magnitude at Bjorken x ~ 0.01. Specifically, the data support 
the hypothesis of a symmetric composition of the light quark sea at low x. The ratio of the strange- 
to-down sea quark distributions is determined to be I.OOIq Jg at absolute four-momentum transfer 
squared Q^ = i.g QeV^ and x = 0.023. 
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A QCD analysis is reported of ATLAS data on inclusive W and Z boson production in pp collisions 
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sensitivity to the light quark sea composition and magnitude at Bjorken x ~ 0.01. Specifically, 
the data support the hypothesis of a symmetric composition of the light quark sea at low x. The 
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Little is known about the strange quark distribution 
in the proton. Flavor SU{3) symmetry suggests that 
the three light sea quark distributions are equal. How- 
ever, the strange quarks may be suppressed due to their 
larger mass. The nucleon strange density plays an im- 
portant role for a number of physics processes, ranging 
from measurements at proton-proton colliders of W bo- 
son production associated with charm jets [I] and of the 
W boson mass P] , to the formation of strange matter [3] 
and neutrino interactions at ultrahigh energies i4j. 

Knowledge of the parton distribution functions (PDFs) 
of the proton comes mainly from deep-inelastic lepton 
proton scattering experiments covering a broad range of 
Q^, the absolute four-momentum transfer squared, and 
of Bjorken x. The PDFs arc determined from data using 
perturbative Quantum Chromodynamics (pQCD). The 
region x < 0.01 is primarily constrained by the precise 
measurement of the proton structure function F2{x, Q^) 
at HERA [5J , which determines a specific combination of 
light quark and anti-quark distributions. However, the 
flavor composition of the total light sea, xS = 2x(u + d + 
s), has not been determined at these x values. 

The strange quark distribution has been accessed in 
charged current neutrino scattering through the subpro- 
cesses W'^s — >■ c and W~s — >■ c. This measurement has 
been made by the NuTeV [5] and CCFR [7] experiments, 
providing information on the strange and anti-strange 
density at x ~ 0.1 and Q^ ^ 10 GeV^. However, the 
interpretation of these data is sensitive to uncertainties 
from charm fragmentation and nuclear corrections. The 
analyses of Refs. [8- 11 suggest strangeness suppression, 
with s/d < 0.5, whereas the analysis of Ref. [12 is con- 
sistent with s/d ~ 1 (unsuppressed strangeness). Recent 
HERMES Kaon multiplicity data [13j point to a strong 
X dependence of the strange quark density and a rather 
large value of x{s -I- s) at a; ~ 0.04 and Q^ ~ 1.3 GeV^. 
However, the interpretation of these data depends on the 
knowledge of the fragmentation of strange quarks to K 
mesons at low Q^. 

In the present letter it is shown that the differential 
measurements of the inclusive W^ and Z boson cross 
sections at the LHC, recently performed by the ATLAS 



collaboration using 35 pb^^ oi pp collision data recorded 
in 2010 [14] , provide new constraints on the strange quark 
distribution at high scale, Q^ ~ M^^, which imply con- 
straints at low Q^ through pQCD evolution. Because 
of the weak couplings of the quarks involved, comple- 
mentary information to F2 is provided which also con- 
strains xYi. A quantity of special interest is the ratio of 
the {W'^ + W^) and Z cross sections which is sensitive 
to the flavor composition of the quark sea [TH [TSl [T5] 
and is rather insensitive to the influence of higher or- 
der pQCD corrections [17 . The inclusive electromag- 
netic and weak DrcU-Yan scattering process is theoreti- 
cally well understood _151 [12] . The parton distribution 
analysis is performed here in next-to-next-to leading or- 
der (NNLO) QCD using the ATLAS data jointly with 
inclusive deep inelastic scattering data from HERA. 

The combined e^p cross section measurements of HI 
and ZEUS [5] cover a kinematic range of Q^ from near 
1 GeV^ to above 10^ GeV^ and of x from ^ 0.6 down 
to lO""*. The ATLAS data access a kinematic range pre- 
scribed by the boson masses, Myi/Zi a-nd the proton beam 
energy, Ep = 3.5 TeV, corresponding to Q^ ~ ^'^w z ^^'^ 
an X range 0.001 < a; < 0.1, with a mean x — AIz/2Ep — 
0.013 for the Z boson. The W^ and Z boson differential 
cross sections have been measured jH] as functions of 
the W decay lepton (e,/i) pseudorapidity, iji, and of the 
Z boson rapidity, yz, respectively, with an experimental 
precision of typically (1 — 2) % in each bin. The abso- 
lute normalization of the three cross sections is known to 
within 3.4 %. Many systematic uncertainties on the mea- 
surements of the W^ and Z boson cross sections are fully 
correlated. These correlations are taken into account in 
the analysis. 

The present QCD analysis uses the HERAFitter frame- 
work [51 [201 m]. The light quark coefhcient func- 
tions are calculated to NNLO as implemented in QCD- 
NUM [22] . The contributions of heavy quarks are calcu- 
lated in the general-mass variable-flavor-number scheme 
of Refs. [23l [24]. The electroweak parameters and cor- 
rections relevant for the W and Z boson production pro- 
cesses are determined following the procedure described 



in Ref. [13], and the results are cross-checked between 
the FEWZ [19] and the DYNNLO [18] programs. The 
HERAFitter package uses the APPLGRID code [25] in- 
terfaced to the MCFM program [5B] for fast calculation 
of the differential W and Z boson cross sections at NLO 
and a iiT-factor technique to correct from NLO to NNLO 
predictions. The data are compared to the theory using 
the x^ function defined in Refs. [27H29]. 

The evolution equations yield the PDFs at any value 
of Q^ given that they are parametrized as functions of x 
at an initial scale Qq. In the present analysis, this scale 
is chosen to be Qq = 1.9 GeV^ such that it is below the 
charm mass threshold rn^. The heavy quark masses are 
chosen to be nic = 1.4 GeV and mi, = 4.75 GeV. The 
strong coupling constant is fixed to as{Mz) = 0.1176, as 
in Ref. [5]. A minimum Q^ cut of Q^i„ > 7.5 GeV^ is 
imposed on the HERA data. 

The quark distributions at the initial scale are repre- 
sented by the generic form 



xq,{x) = A,x'''il-x)'^^P,ix) 



(1) 



where Pi{x) denotes polynomials in powers of x. The 
parametrized quark distributions, xqi, are chosen to be 
the valence quark distributions (xuy, xdy) and the light 
anti-quark distributions (xu, xd, xs). The gluon dis- 
tribution is parametrized with the more flexible form 
xg{x) = AgX^!>{l - x)^i'Pg{x) - A'gX^'!i{l - x)'^'s , where 
C'g is set to 25 to suppress negative contributions at high 
X. The parameters Au„ and Ad^ are fixed using the quark 
counting rule and Ag using the momentum sum rule. The 
normalization and slope parameters, A and -B, of u and 
d are set equal such that xu = xd at x ^ 0. Terms are 
added in the polynomial expansion Pi{x) only if required 
by the data, following the procedure described in Ref. [5] . 
This leads to one additional term, Pu^{x) = 1-1- E^^^x'^ . 

Two types of NNLO fit, termed epWZ, are performed 
with different treatments of strangeness. First, the 
strange quark distribution is fully coupled to the down 
sea quark distribution and suppressed by fixing sjd = 0.5 
at the initial scale Qg ("fixed s fit") as suggested by 
Refs. [3 [8l-[TT]. In a second fit, xs is parametrized as 
in Eq.[Tl with Pg = 1 and Bg = B^, leaving two free 
strangeness parameters, Ag and Cs ("free s fit"). By 
default it is assumed that xs = xs. 

Both fits result in good overall x^/^df values of 
546.1/567 with 13 free parameters, for fixed s, and of 
538.4/565 with 15 free parameters, for free s. For the 
fixed s fit, the partial x^ of the ATLAS data is 44.5 
for 30 data points. This improves significantly to 33.9 
for the fit with free s. This fit determines the value of 
Ts — 0.5(s + s)/d to be 

rs = 1.00±0.20cxp±0.07modl°;i°part°°?as±0.08th, (2) 

at Qo and x = 0.023, the x value, which corresponds to 
X = 0.013 at Q^ = M| as a result of QCD evolution. 



The combined result is r. 



i.oo: 



-0.25 
-0.28- 



The uncertainty of rg (Eq.[2]) is dominated by the ex- 
perimental (exp) uncertainty, which is mostly driven by 
the statistical and systematic uncertainties of the W and 
Z cross section measurements. The model (mod) uncer- 
tainty includes effects due to variations (1.25 < rric < 
1.55 GeV and 4.5 < nib < 5.0 GeV) of the charm and 
beauty quark masses following Ref. |30| , of the minimum 
Q'^ cut value (5 < Q^m < 10 GeV^), and the value of the 
starting scale (Qg lowered to 1.5 GeV'^). The largest con- 
tribution to the model uncertainty of ±0.05 comes from 
the variation of the charm quark mass. The parametriza- 
tion (par) uncertainty corresponds to the envelope of the 
results obtained with the polynomials Pi, in Eq.lll ex- 
tended by one or two terms, resulting in somewhat differ- 
ent parton distributions with similar x^ as for the nom- 
inal fit. The parametrization uncertainty also includes 
a fit with Bg free. The a^ uncertainty corresponds to 
a variation of as{Mz) from 0.114 to 0.121. Finally, a 
theoretical (th) uncertainty is assessed by comparing the 
DYNNLO and FEWZ predictions on the Z, W+ and W' 
fiducial cross sections, which agree at the level of 0.2, 0.5 
and 1.0%, respectively. In addition, remaining missing 
pure electroweak corrections may alter the QCD predic- 
tions at the per mille level. Both effects are well covered 
by an uncertainty of 1% on the W/Z cross section ratio 
and this results in a theoretical uncertainty on r^ of 0.08. 

The fits impose small shifts, typically much smaller 
than one standard deviation, on the correlated system- 
atic uncertainties of the data. The global normalization 
is observed to be shifted upwards for both fits by about 
the size of the luminosity measurement uncertainty. The 
W^ and Z cross section measurements are compared in 
Fig. [T] to the NNLO fit results, after these shifts are ap- 
plied to the predictions. Also shown are the ratios of the 
fits with free s and with fixed s. It is apparent that the 
enhanced strange quark fraction in the free s fit has no 
significant effect on the prediction of the rji distributions 
for both the W^ and W~ decay leptons, while it leads 
to an improvement in the prediction of the yz distribu- 
tion. An improvement is also observed in the description 
of the ratio of the {W~^ + W~) to the Z boson cross sec- 
tions in the fiducial phase space. This is predicted to be 
11.10 in the fit with fixed s, while the measured value of 
10.70 ± 0.15 is almost exactly reproduced in the fit with 
free s, which gives a value of 10.74. 

In order to check the robustness of the present result 
for Ts , a series of cross checks is performed. A fit without 
allowing an adjustment of the correlated errors yields a 
value of Tj = 0.97±0.26oxp, in good agreement with Eq.[2] 
A fit with identical input parameters is repeated at NLO 
and also yields a consistent result: r^ — 1.03 ± 0.19cxp. 
If this NLO fit is performed with a massless heavy quark 
treatment then r^ = 1.05 ± 0.19exp is obtained. In a 
separate NLO study, the constraint {xu — xd) — >■ for 
X — >■ is relaxed. The xd{x) distribution is found to 
be consistent with xu{x), albeit with large uncertainties 
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FIG. 1. Differential da/d\rii+\ (left) and da/d\rji-\ (middle) cross section measurements for W -^ iv and dcr/djj/z| cross 
section measurement for Z -^ (.1 (right). The error bars represent the statistical and uncorrelated systematic uncertainties 
added in quadrature while the theoretical curves are adjusted to the correlated error shifts (see text). The NNLO fit results 
with free and fixed strangeness are also indicated, and their ratios are shown below the cross section plots. 



(~ 15% at a; ~ 0.01 and Qq). The fraction of strangeness 
is again consistent with unity, r^ = 0.96±0.25oxp. FinaUy 
the data are fitted, to NNLO, with separate strange and 
anti-strange normahzations. The resulting value of r^ is 
consistent with unity and the ratio s/ s is 0.93 ± 0.15exp 
at a; = 0.023 and g^ = Ql. 

W, Z cross section measurements performed at the 
Tevatron may potentially have sensitivity to rg similar 
to that of the ATLAS data. A NLO fit to the HERA 
with the CDF W asymmetry [31] and Z rapidity [32] 
data gives r^ = 0.66 ± 0.29oxp at a mean x of about 
0.081. This is consistent within uncertainties with both 
suppressed strangeness and with the present result. A 
NLO fit to the combined HERA, ATLAS and CDF data 
yields r^ = 0.95 ± 0.17oxp. 

The provision of the full differential cross sections for 
both W~^, W~ and Z boson production, besides the ep 
cross sections, is essential for the determination of xs: if 
the ATLAS Z cross section data are fitted together with 
the ATLAS W charge asymmetry data, rather than with 
the separate W^ and W~ cross section measurements, a 
less precise result is obtained with r,, = 0.92 ± 0.31cxp. 

In Fig.[2]the present result for r^ is compared with pre- 
dictions obtained from four global PDF determinations. 
The CTIO (NLO) [H] determination gives a large frac- 
tion consistent with the present result. On the other 
hand, the MSTW08 [8] and ABKM09 ^9 determina- 
tions give a much lower value of r^ ~ 0.5, and the 
NNPDF2.1 [TUldl] result of r^ ~ 0.25 is even lower. 

The enlarged fraction of the strange quark sea leads to 
a decrease of the down and up quark sea densities at the 
initial scale Qq, because xs, xd and xu are tied together 
at low X by the precise F2 data. In compensation for the 
increase of xs, the xd and xu distributions are dimin- 
ished by ~ 10%. The total sea, a;S, is correspondingly 
enhanced by ~ 8%, as illustrated in Fig.jsj 
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FIG. 2. Predictions for the ratio r^ = 0.5(s -f s)/d, at 
Q2 = 1.9GeV^ X = 0.023. Points: global fit results us- 
ing the PDF uncertainties as quoted; bands: this analysis; 
inner band, experimental uncertainty; outer band, total un- 
certainty. 



The result on r^, Eq.[2j evolves to 

Ts = 1.00±0.07oxp±0.03mod+Jj°^par±0.02Qs±0.03tti (3) 

at Q^ = Af| and x = 0.013, corresponding to a value 



of r, (0.013, Af|) = 1.00 



+0.09 
-0.10' 



which is more than twice 



Ml 



as precise as at the initial scale Qq. Uncertainties are 
smaller at Q^ = M^ because the gluon splitting proba- 
bility into qq pairs is flavor independent, thus reducing 
any initial flavor asymmetries. This also causes rs to in 
crease from 0.5 at Qq to a value of about 0.8 at Q^ 
in the fixed s fit. 

In summary, a NNLO pQCD analysis is performed of 
the first differential ATLAS W^ , Z pp cross sections with 
HERA e^p data. The W, Z measurements introduce a 
novel sensitivity to the strange quark density at a; ^ 0.01, 
which is exploited here for the first time. The ratio of 
the strange to the down sea quark density is found to be 
Vg = l.OO^Q 28' ^^ Bjorken x = 0.023 and the initial scale 
of the QCD fit Ql = 1.9 GeV^. This is consistent with 
the prediction that the light quark sea at low x is flavor 
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FIG. 3. Distribution of the light sea quarks, a;S = 2x{u + 
d + s), in the NNLO analysis of HERA and ATLAS data with 
a fixed fraction of strangeness (lower, green curve) and with a 
fitted fraction of about unity (upper blue curve). The bands 
represent the experimental uncertainties. 



symmetric. A consequence of this initial observation is 
that the total sea, xT, = 2x{u + d + s), is enhanced by 
about 8%, as compared to the result when the strange 
quark is suppressed to half the magnitude of the down 
sea quark. 
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